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ABSTRACT 

We investigate the relationship between black hole mass and bulge luminosity for AGNs with reverberation- 
based black hole mass measurements and bulge luminosities from two-dimensional decompositions of Hubble 
Space Telescope host galaxy images. We find that the slope of the relationship for AGNs is 0.76-0.85 with 
an uncertainty of ^0.1, somewhat shallower than the Mbh oc i relationship that has been fit to nearby 

quiescent galaxies with dynamical black hole mass measurements. This is somewhat perplexing, as the AGN 
black hole masses include an overall scaling factor that brings the AGN Mbh - o"* relationship into agreement 
with that of quiescent galaxies. We discuss biases that may be inherent to the AGN and quiescent galaxy 
samples and could cause the apparent inconsistency in the forms of their Mbh -i^buige relationships. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: photometry — galaxies: Seyfert 



1. INTRODUCTION 

Most galactic bulges are now believed to harbor a massive 
black hole. For AGNs, the evidence of the black hole is obvi- 
ous from the activity. However, even in quiescent galaxies, the 
effect of the black hole can be detected from stellar and gas 
kinematics near the nucleus. What remains to be determined 
is the formation mechanism for these massive black holes and 
their role in shaping and responding to the evolution of their 
host galaxies. 

In an early revie w on the subject of q u iescen t galaxy 
black hole masses, iKormendv & Richstond (1 19951) pointed 
out that the estimated central black hole masses of eight 
galaxies seemed to show a correlation with the host galaxy 
bulge luminosit i es (or, equivalently, bulge stellar masses). 
iMagorrian et al.l (Il998h later investigated a sample of 32 
nearby galaxies and their black hole masses, and confirmed 
that the estimated black hole mass in each galaxy was indeed 
proportional to the luminosity (or mass) of the host galaxy 
bulge, albeit with a scatter about the fit of approximately 
±0.5 dex. Later studies claimed, in some cases, that there 
may be a difference in the black hole-bulge relationship for 
quiescent and active galaxie s, or even for Seyfert galaxies and 
quasars (e.g.. I Wandelll 1 999l) . Subsequent investigations seem 
to have decreased these discrepancies through more sophisti- 
cated techniques of measuring the bulge lu minosity such as 2- 
dimensional image decompositions (e.g., iMcLure & DunlopI 
l200lHWandelll2002l rhe reafter W021) or dyn amical modehng 
of the host galaxy (e.g.j Haring & Rixll2004l) . 

We have recently completed a Hubble Space Telescope 
campaign to image the host galaxies of AGNs with black 
hole masses from reverberation-mapping. The images were 
acquired for the purpose of decomposing the surface bright- 
ness profiles of the host galaxies and creating "nucleus- 
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free" images from to measure the host-galaxy starlight con- 
tributions to ground-based spectroscopic lum inosity measure- 
ments of the AGNs (results described by B entz et alj|2006l 
l2008l) . Combining the bulge luminosities estimated from the 
surface brightness decompositions of these high-resolution 
images with the recently updated and homogeneously an- 
alyzed database of reverberation masses for these objects 
(iPeterson et al.ll2004t iGrier et al.l l2008i) . we re-examine the 
relationship between black hole mass and host galaxy bulge 
luminosity for nearby AGNs. Throughout this work, we will 
assume a standard flat ACDM cosmology with JIb = 0.04, 
f^DM = 0.26, f^A = 0.70, and //q = 70kms-' Mpc"'. 

2. BULGE LUMINOSITIES AND BLACK HOLE MASSES 

We investigate here the sample of AGNs with black hole 
masses measured using the variability technique known as re- 
verberation mapping. The details of the host galaxy imag- 
ing and sur face brightness decompositions are described by 
iBentz et aP (2006, 2008), but we include a short summary of 
the relevant details here. The majority of the objects in the re- 
verberation sample were imaged with the Advanced Camera 
for Surveys (ACS) High Resolution Channel (HRC) through 
the F550M filter Unfortunately, ACS ceased function- 
ing before the observations were completed, so five objects 
(PG0026, PG1307, PG 1426, PG1617, and Mrk509) were 
imaged with the Wide Field Planetary Camera 2 (WFPC2) 
through the F547M filter Exposure times for all observa- 
tions were graduated so that the saturated pixels in the nu- 
cleus region of the long exposures could be corrected using 
an unsaturated image. After correcting saturated pixels and 
cleaning cosmic rays, individual images were co-added and 
corrected for distortion if necessary. The surface brightness 
profiles of the AGN host galaxies were modeled with the two- 
dimensional image decomposition program Galfit (Peng et aO 
l2002h . To avoid overcorrecting the AGN luminosities for host 
galaxy starlight, we made conservative host galaxy model fits 
which may somewhat underestimate the true brightness of the 
various galaxy components. 

For most of the reverberation-mapped objects, the HST im- 
ages are sufficient for a full decomposition. However, for the 
nearby and spatially-extended NGC objects in the sample, the 
HST images do not provide the necessary wide field of view 
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Note. — For those objects with two bulge luminosities listed, the first is the luminosity of the large st-scale non-di s k com ponent, and the second is 
the luminosity of all non-disk components, including bars or inner bulges. Black hole masses are from lPetersonetalJ <2004l) , except for PG2130+099 
IGrier et al.ll2008l) . 



(FOV) to measure the sky background, which affects the accu- 
racy to which we can constrain the bulge and disk parameters. 
Our current surface brightness models for the HST images of 
the NGC galaxies include components that cannot be resolved 
in the ground-based images due to their small effective radii 
(<3.5" in all cases). The typical seeing in the ground-based 
images of ^ 2" blurs the central PSF together with these com- 
pact components and with the bulge itself, which has an effec- 
tive radius of only ^ 10" for these objects. The difficulties in 
disentangling the nuclear structure in the ground-based im- 
ages and the fact that the distortion-corrected HRC FOV is, at 
best, only twice the effective radius of the bulge and further- 
more doesn't allow measurement of the sky, creates enough 
uncertainty in the bulge parameters for the NGC objects with 
reverberation masses that we exclude them from further anal- 
ysis at this time. Additional images with an intermediate FOV 
are necessary to fully constrain the decompositions of these 
galaxies. We also exclude IC4329A because of its uncertain 
black hole mass and because it is a dusty edge-on galaxy with 
an unreliable surface brightness decomposition. Table 1 lists 
the 26 objects examined by this paper. 

The paramet ers of the surface b rightness decompositions 
are included in Bentz et al] (|2008|) . We take the component 
with the largest effective radius (other than any exponen- 
tial disk component) as the "bulge" for each of these galax- 
ies, except for Mrk79 where there is clearly an extended 
bar-like component. For ellipticals, the bulge magnitude is 
th e total galaxy magni tude. We followed the prescriptions 
of ISirianni et al.l (l2005l) for converting from the space tele- 
scope AB A magnitude system to Galactic extinction-corrected 
magnitudes in the Vega system, which we list in Table 1 as 
niHST- The V -HST color for each object was c alculated with 
synphot in IRAF using the bulge template of iKinnev et al.l 



(Il996l) . We include in Table 1 the extinction-corrected in- 
tegrated bulge magnitude in the Vega system for each AGN 
host galaxy, the corresponding apparent and absolute John- 
son V magnitudes, and the V-band luminosity from the stan- 
dard relation logLy/i^o = 0.4(-Mv + 4.83). We assume an 
uncertainty of ±0.2 dex for the bulge luminosities. Four of 
the galaxies had an additional surface brightness component 
besides a bulge or exponential disk. Imaging alone is not suf- 
ficient to ascertain whether these components contribute to the 
total bulge brightness or are dynamically distinct, so we also 
list in Table 1 the bulge luminosity for those objects including 
the contribution from the additional component. 

Figure 1 compares the (single component) bulge luminosi- 
ties and black hole masses for the sample of objects that are 
common to both this work (filled circles) and the Wandel 
(I2OO2I) study (open circles). There is a clear trend between 
black hole mass and bulge luminosity for th e objects in this 
study, whereas the bulge luminosities from IWandell (l2002h 
have a range of ~ 2 dex in black hole mass over a Um- 
ited ^ 1.5 dex range in luminosity. It is perhaps unsurpris- 
ing that a correlation is more apparent now, as the study by 
Wandel (2002) is a compilation of inhomogeneous data from 
the literature. The bulge luminosities were primarily from 
iMcLure & DunlopI (^001) and include measurements from 
ground-based photographic plates and CCD imaging, as well 
as saturated and unsaturated //5r imaging, all in various pass- 
bands. Many of the studies included in the compilation did not 
use the same cosmologies and distances for the objects or the 
same analysis techniques in determining the black hole mass. 
All of the reverberation-based black hole masses have since 
been h omogeneously analyzed and updated by Peters on et alJ 
(12004 . and PG2130H-099 which has an updated mass from 
the analysis of a new reverberation data set by iGrier et al] 
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Fig. 1. — Compaiison of black hole mass and V-band bulge luminosity 
values from Wandel ( 2002) (open points) and this work (filled points) for 
objects that are common to both. The new values cover a range of 2.5 dex in 
luminosity and show a clear trend, while the values from Wandel ( 2002) have 
a range of ~ 2 dex in mass within a limited ~ 1 .5 dex range in luminosity. 

diOOS). Fi sure 2 shows the black hole masses versus the host 
galaxy V-band luminosities for the full sample of 26 objects 
included in this study. 

3. THE BLACK HOLE MASS - BULGE LUMINOSITY 
RELATIONSHIP 

We employed two independent fitting routines i n our exam- 
inatio n of the MfiH-ibuige relationship: FITEXY dPress et alj 
|1992|) . which estimates the parameters of a straight-line fit 
throug h the data including error s in both coordinates; and 
BCES jAkritas & Bershadvll 19961) . which accounts for the ef- 
fects of errors on both coordinates in the fit using bivariate 
correlated errors and a component of intrinsic scatter. FI- 
TEXY numerically solves for the minimum orthogonal us- 
ing an interative root-finding algorithm and is a "symmetric" 
algorithm in that it does n ot assume a dependent and indepen- 
dent variable. Following iTremaine et alj (l2002h . we include 
an estimate of the fractional scatter, in this case the fraction of 
the Mbh measurement value (not the error value) that is added 
in quadrature to the error value to obtain a reduced of 1.0. 
While BCES also accounts for intrinsic scatter, it does not 
provide a measure of it. We adopt the bootstrap of the BCES 
bisector value with = 1000 iterations. Fits of the form 
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were performed utilizing both the single-component bulge 
luminosities and the multiple-component bulge luminosities 
and are presented in Table 2. The powerlaw slope ranges from 
0.76 ± 0.08 to 0.85 ± 0. 11 depending on the definition of the 
bulge luminosity and the specific fitting routine utilized. We 
take the BCES fit nearest the middle of this range, with a slope 
of 0.80 ± 0.09, as the "best" fit. 

For comparison, we fit the quiescent galaxy Mbh ~ i'buige 
relationship using the sample of near by galaxies with dynam- 
ical black hole mass measurements (iFerrarese & FordI 120051: 



9 - 



cn 
O 




This Worl< 
Wandel 2002 

, I 



1 1 



12 



log U (Lg) 



Fig. 2. — The MflH "^buloc relationship for AGNs with reverberation-based 
masses and bulge luminosities from two-dimensional decompositions of HST 
host-galaxy images. The solid line is the BCES fit w ith a slope of q = 0.80 ± 
0.09. The dashed line is the fit from lWandell 120021) to his sample of broad- 
line AGNs and has a slope of o = 0.90 ± 0. 1 1 . 

FF05). We restricted the sample to ellipticals, both to cir- 
cumvent the need for bulge-disk decompositions and be- 
cause ellipticals are reporte d to show less scatter ab out the 
MBH-ibuige relationship (cf. iMcLure & Dunlopl200lh . Bulge 
magnitudes were converted to V-band using a typical ellipti- 
cal galaxy color of B-V = 0.9. The fitting results are pre- 
sented in Table 2. We also fit the quiescent galaxy relation- 
ship excluding CygnusA and NGC5845, both of which are 
known to deviate significantly (FF05). 

Figure 3 shows the Mbh - ^buige relationship for the FF05 
ellipticals compared to the "best" fit for the AGNs in this 
study. The slope of the quiescent galaxy relationship is 
steeper, although the severity of the discrepancy depends on 
the specifics of the fitting routine used, as the best-fit slope 
from BCES is a = 1.43 ±0.21 versus a = 1.11 ±0.21 from 
FITEXY. The fits presented here for the quiescent ellipti- 
cal galaxies are slightly steeper than that reported by FF05 
(a = 1.05 ±0.21) and are in reasonable agreement with the 
slope of 1.40 ±0.1 7 for the quiescent galaxy fit reported by 
iLauer e t aT ('2007a'). Both of these studies included lenticular 
and spiral galaxies with dynamical masses in their fitting sam- 
ples. ILauer et al.l (I2007al) compiled decompositions from the 
literature, while FF05 assumed specific ratios of Lbuige/^-totai 
based on morphological type.^ 

4. DISCUSSION 

The compendium of AGN black hole masses and bulge 
luminosities presented here offers several advantages over 
previous measurements compiled from the literature (i.e., 
'Wanderi999l iMcLure & Dunlop"2001). The masses 

result from a homogeneou s analysis of reverberation-mapping 
data jPeterson et alj|2004h . in contrast to many recent stud- 

^ The relative contribution of the bulge is known to va ry substantially 
within a single morphological type. Figure 6 of lKenll 119851) shows, for ex- 
ample, that B/r ranges from 0.3 — 1.0 for SO galaxies). 
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TABLE 2 

Fits to the Mbh -^buise Relationship 
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Sample 


K 


a 


Scatter" 


BCES 


AGNs 

AGNs (+ extra components ) 

FF05 Ellipticals 

FF05 Ellipticals (— outliers) 


-0.02 ±0.06 
-0.07 ± 0.08 
0.42 ±0.12 
0.30±0.10 


0.80 ±0.09 
0.85 ±0.11 
1.43±0.21 
1.42 ±0.24 




FITEXY 


AGNs 

AGNs (+ extra components) 

FF05 Ellipticals 

FF05 Ellipticals (— outliers) 


-0.05 ±0.06 
-0.13 ±0.06 
0.15±0.11 
0.14±0.11 


0.76 ±0.08 
0.80 ±0.09 
1.11±0.21 
1.18±0.19 


0.38 
0.44 
0.73 
0.64 



" The fractional scatter, quantified as the fraction of the measurement 
value of Mbh that must be added in quadrature to the eiTor value in 
order to obtain a reduced x' of 1.0. 

ies of AGN black hole - bulge relationships where black 
hole masses are inferred from single-epoch s pectral mea- 
surements (e.g., [Oriiniitin llOoCllImiiti^ The 
bulge luminosities in this work are estimated from two- 
dimensional surface brightness decompositions of unsatu- 
rated high-resolution space-based images taken with the same 
instrument and the same filter The exceptions are the five 
objects that were imaged through the F547M filter using 
WFPC2, but the bandpass is very similar to the ACS F550M 
filter employed for the other objects (Ac(F550M) = 5580 A 
versus Ac(F547M) = 5483 A, and AA(F550M) = 547 A versus 
AA(F547M) = 483 A). 

Unfortunately, there is not a similarly consistent sample 
of high-quality images from which the bulge luminosities 
can be estimated for the quiescent galaxies with dynami- 
cal black hole masses. While high-quality observations do 
exist for these nearby and well-studied galaxies, they have 
no t been obtained in a uniform fashion. The recent work 
by iGrahamI (l2007h attempts to compensate for the different 
methods and ana lysis techniques employed in several pub- 
lications dMcLu re & Dunlop 2002 w ith an updated cosmol- 
ogy presented in .McLure & DunlopI l2004t iMarconi & Hunt! 
120031: lErwin et alJl2004 . all of which arrive at different val- 



ues for the slope of the quiescent galaxy Mbh - ^buige relation- 
ship and/or different black hole masses predicted for a specific 
bulge luminosity. iGrahamI carefully updated and revised the 
samples of objects included in these studies and found that 
the differences of the best-fit parameters found by each study 
are mitigated, the scatter in the measurements is significantly 
decreased, andthatMBH ocL' Interestingly, he finds asome- 
what shallower slope, a w 0.75-0.88, when bulge luminosi- 
ties are estimated by two-dimensional decompositions of B- 
band images and corrected for inclination-dependent dust ex- 
tinction in the host galaxy disks. Only 13 objects are included 
in that par ticular analy sis (an updated form of the study pre- 
sented by Erwin et al .' 2004) but it is intriguing nonetheless in 
its close agreement with the fit that we find for the AGNs. 

While we expect that the AGN bulge luminosities in this 
may be slightly underestimated based on the conservative 
galaxy fits we employed, there is also a small A:-correction 
introduced, as the average redshift of the 26 AGNs is z w 0.1. 
The portion of the SED observed through the medium-band V 
filters employed in the HST observations is fainter 0. 1 dex 
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Fig. 3. — The MBH~^'bulgc relationship for quiescent galaxies as deter- 
mined for the elliptical galaxies from Femrese & Ford (2005; filled trian- 
gles) compared to the Mbh ~ ^-bulgc relationship for AGNs from this work. 
The slope of the relationship appears shallower for AGNs than for quiescent 
galaxies. The solid line is the same fit displayed in Figure 2 for the sample 
of AGNs in this work, and has a slope of a = 0.80 ± 0.09. The short-dashed 
line is the BCES fit to the elliptical galaxies and has a slope of 1.43 ±0.21, 
while the long-dashed line is the FITEXY fit to the elliptical galaxies and has 
a slope of 1.11 ±0.21. 



in luminosity) than if the galaxies were at z = 0, assum- 
ing thei r stellar populations resemble that of the bulge tem- 
plate of iKinnev et alj (11996'). Accounting for any such bi- 
ases in the galaxy fitting or colors would intensify the appar- 
ent differences between the slopes measured for the AGNs 
and quiescent galaxies. The black hole masses for the AGNs 
have already been scaled so that their Mbh - relationship 
is brought into agreement with the quiescent Mbh - cr* re- 
lationship fit to the same sample of quiescent galaxies in- 
cluded in the FF05 s t udy. The differences may be miti- 
gated if iMarconi et al.l (l2008h are correct in their suggestion 
that neglecting radiation pressure leads to systematic under- 
estimation of black hole masses from reverberation mapping 
data (see, however, Netzer 2008). 

A s ep arate concern h as be en raised by Yu & Tremaind 
(120021). iBernardi et al.l (l2007h . iLauer et al.1 



- Y7 r (I2007al) . and 

iTundo et al.l (l2007l) . who present an apparent disagreement 
between the quiescent galaxy Mbh - and Mbh - ^buige re- 
lationships and suggest that the quiescent galaxy sample is 
biased towards galaxies with overly la rge velocity d ispersions 
for their luminosities (see, however, IGrahamI l2008i who ex- 
amines the role of bars in this issue). There is no reason to 
suspect the AGN sample of having the same bias, as the mass 
measurements are made using flux variability techniques and 
not dynamical techniques. Such a bias may help explain why 
some of the quiescent galaxies at the high luminosity end have 
black hole masses that are more than an order-of-magnitude 
larger than the active galaxies, although it may not completely 
resolve this disparity. 

Finally, there may be no reason to expect that the Mbh - 
i'buige relationship is the same for the AGNs and quiescent 
galaxies in these samples, as there is only a modest number of 
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objects in each sample and selection effects lik ely play an im- 
portant role on both sides (iLauer et alJl2007bl) . Clearly, there 
remain several areas that are in need of investigation, any of 
which may shed light on the apparently inconsistent fits to the 
MBH-i-biiige relationship for AGNs and for quiescent galaxies. 
As the Mbh - ^buige relationship is an important and widely 
used means of est imating black hole masses throughout cos - 
mic history (e.g.. iMarconi et alj 120041: IShankar et alj|2004 . 
an accurate characterization of this relationship is necessary 
for understanding black hole growth and evolution as well as 
the interplay between black holes and their host galaxies. 

5. SUMMARY 

We have presented an updated version of the AGN Mbh - 
ibuige relationship using the datab ase of homoge neously an- 
alyzed reve rberation masses from Peterson et al.l (2004) and 
iGrier et al.1 (2008; PG 2 130+099) and the two-dimensional 
surface brig htness decornpositio ns of the AGN host galaxies 
described by iBentz et al.l OOOSh . We find a strong correlation 
about the relationship for the 26 AGNs included here, with a 
best-fit powerlaw slope of 0.80 ±0.09. This is somewhat shal- 
lower than the best-fit slope for quiescent galaxies (a « 1.0), 
even though the AGN black hole masses have been scaled to 
bring the AGN and quiescent galaxy Mbh - cr* relationships 
into agreement. There appear to be many systematics in both 
the AGN and quiescent galaxy samples that must be investi- 
gated in order to more completely understand this important 
relationship. 

Our future plans include investigating the biases in the AGN 
sample and extending the range of the relationship for AGNs. 
We have an //ST Cycle 17 program to image the NGC objects 
that were excluded from this particular work with the Wide 



Field Camera 3 through the F547M filter. These observations 
will provide us with the intermediate FOV images necessary 
for accurate decompositions of those galaxies, enabling us to 
include them at the low-mass end of the MBH-^buige relation- 
ship for AGNs. Recent reverberation-mapping experiments 
that were carried out at MDM Observatory (spring 2007) and 
Lick Observatory (spring 2008) focusing on AGNs with black 
hole masses in the range 1 x 10^-5 x IO^Mq (Denney et al. 
in preparation, Bentz et al. in preparation) show promise in 
further extending the range and coverage of the Mbh - i'buige 
relationship for AGNs at the low-mass end. 
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